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Abstract 

We analize two reactions: the electroproduction at small x and the Drell-Yan pro- 
duction by protons on nuclei. We show that nuclear shadowing in electroproduction is 
larger if the virtuality of intermediate hadrons is properly taken into account. In the 
Drell-Yan production, initial state interactions of fast partons are included. With this 
input, the excess pion contribution improves the description of data for both reactions. 



1. Following the conventional picture of nuclei, the meson exchange currents result in a 
certain amount of excess pions in nuclei. This fundamental nuclear property was questioned 
by some recent experiments. The behaviour of nuclear structure functions at low and 
a lack of nuclear cross section enhancement for Drell-Yan (DY) production by protons |2|] 
have led several authors (see, e.g., the Ref.0) to the conclusion that there are no excess 
pions in nuclei. The present analysis of these reactions is based on some new conjectures. In 
electroproduction, it is shown that virtual hadrons interact only inelastically with targets. 
As the result, the nuclear shadowing can be significantly stronger than it has been believed 
so far. In DY production by protons, as well as in other hadroproduction reactions Q, initial 
state interactions of fast partons are taken into account. With this input, the data for both 
reactions are much better described with the excess pion contribution included. 

2. In the impulse approximation, the quark (antiquark) distribution in nuclei can be 
represented as a convolution of nucleon and excess pion quark distributions, 

J v v j v y 

where f§(y) and f A (y) are the nucleon and excess pion distribution functions in nuclei. 
The flavour indices and the Q 2 — dependence are implied here and hereafter. Note that 
we denote by f A {y) the renormalized excess pion distribution which preserves the total 
probability conservation ||. Therefore the functions f3(y) and f^(y) are not independent 
but are constrained by the momentum sum rule 

/ dyyf£(y) + J dyyf^y) =< y A > + < y* >= 1. (2) 

One way to find q A (x) is to calculate first f A {y) using the elementary nN interaction and then 
to determine f^(y) from (2), adopting some model for the nucleon momentum distribution. 
In earlier works of this kind|| the unrenormalized pion distribution was used for f A (y), 
what led to an overestimation of the pionic contribution to q A (x). This was corrected in 
more recent papers^, but there is a significant model dependence of calculations, which 
start from the elementary rrN interaction. 

Another approach was suggested in Ref.||. Since f^{y) is an invariant function (y is the 
ratio of two scalar products), it can be calculated in any frame. In the nucleus rest frame, 
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/jv(y) is determined by the energy and momentum distribution of nucleons, which can be 
measured in independent experiments]^]. In this case the magnitude of the EMC-effect is 
proportional to (1— < >)/mjy = e/^iv, where e is the average nucleon separation energy. 
For heavy nuclei the data give e >40 MeV. The A-dependence of e was discussed in Ref . |TT| . 
When f^(y) is determined, the first moment of the excess pion distribution is fixed by (2) 
and a detailed form of f^(y) can be found by using some realistic form of pion momentum 
distribution, e.g. the one from Ref.fll^]. Here we calculate f^(y) with this method. 

3. It was shown|8[] that the main part of the EMC-effect can be explained by the binding 
correction of nuclear structure functions. We here focus on the low-x region and show that 
the nuclear shadowing and antishadowing can be explained provided (a) the virtuality of 
intermediate mesons is taken into account while considering their multiple interactions in 
nuclei and (b) the contribution of excess pions is added. We use the generalized vector 
meson dominance (GVMD) model, but the above conclusion is valid for some other hadron 
dominance models, adopting the Glauber theory |T0| to describe multiple interactions of 
virtual hadrons. In the GVMD model, the structure function is given by (see, e.g., the Ref. 



0) 



t-iGV M D / \ Q 2 f°° i 2 /^n^ 2 ) t 2 \ fo\ 

F 2 ( X ) = / d/1 - —r= (Jtot{V ,S), (3) 
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where II(/i 2 ) is the spectral function of mesonic states with the mass /i, a to t is the total 
cross section of meson interaction with a target. Total cross sections are connected to 
forward scattering amplitudes by the optical theorem, a to t = I m [f(s,t = 0)]/s (we use the 
normalization < p\p' >= (2Tr) 3 2p 5(p — p')). In some other models [fL4ll , the parton-target 
amplitudes have been used instead. In the independent particle model, the Glauber result, 
adopted by many authors, is given by 

f A (s, t) = -2si ± d J db C fN{ if A T{b) r, (4) 

where T(b) is the optical thickness of the nucleus and b is the impact parameter. For 
Im[f N ]/Re[f N ] ^> 1 and A 1, (4) can be rewritten as 

af ot = 2 fdb(l-e-<^). (5) 



We now show that (4) and (5) cannot be used for virtual hadrons participating in electro- 
production. 

In Glauber theory it is assumed that for fast projectiles the overall modification of the 



wave inside the nucleus is the factor (JTOJ, [HJ 



l-T A (b) = U[l-T N (b-s i )}, (6) 

i=i 

where T^lb — s$) is the profile function of a nucleon with the transverse coordinate S{. 
This equation is also valid for virtual projectiles if their propagation length A satisfies the 
condition A ^> r A , where r A is the nucleus size. A correction of (6) in the case A ~ r A has been 



considered in Ref.[16 and will be taken into account below. The equation (4) follows from 



(6) if the elastic scattering amplitudes are expressed in terms of profile functions as[|T(], [OJ 

f NA (s,t) = ^ [ d 2 be-™T N!A (b), (7) 
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where k is the projectile three-momentum in the lab. frame. The last equation describes the 



Fraunhofer diffraction of waves [17 1, if the diffraction is observed at large distances z from 



the target. The condition for the Fraunhofer diffraction (and for the validity of (7)) is [17 



z > L, L = 2krf, 



where L is the healing length of shadow [TjJ (L characterizes the duration of the elastic 
scattering process) and r t is the target transverse size. At z ~ L the Frenel diffraction, 
which is weaker than the Fraunhover diffraction, takes place instead and at z <C L there is 
no diffraction |I7|. But for virtual projectiles the accessible longitudinal distances are only 



of the order of A. Therefore the standard derivation of (4) doesn't hold if |18 



A < L. (9) 

In that case, a projectile vanishes before the diffraction takes place. Let us obtain the same 
result using a more formal language. 

For stable projectiles, the Glauber result for elastic and inelastic cross sections reads 

a el = f d 2 b\T(b)\ 2 , a inel = f d 2 b{\ - |1 - r(6)| 2 ). (10) 



As it follows from this equation, both a d and a ine i are non- vanishing if F(b) 7^ 0. However, 
this is not the case for virtual projectiles and we now show that a e i vanishes under the 
condition (9), but Oimi does not. The total interaction probability is given by the imaginary 
part of the forward projectile-target amplitude f(s, t = 0) =< k, P\T\k, P > where k, P are 
projectile and target momenta and the operator T describes all projectile-target interactions. 
The imaginary part of the amplitude can be represented as the sum of elastic and inelastic 
channels, 

Im < k, P\f\k, P >= ^ Yl 5A ( k n + Pn~k-P)<k, P\f\k n , P n >< k n , P n \f\k, P > 

+ E^ 4 ( p * -k-p)<k, p\f*\i >< i\f\k, p >, (ii) 

I 

where \k n , P n > are all possible states of initial particles and \i > are all other states, 
\i >t^ \k n ,P n >. The intermediate states \k n ,P n > and \i > contribute to Im[f(s,t = 0)] 
if they contain only free particles. In fact, to take the imaginary part of a diagram is 
equivalent to cut this diagram and to put the cut lines on the mass shell. In the case of 
virtual hadron-nucleus scattering, P 2 = m\ but k 2 ^ fi 2 , where \i is the hadron mass. It 
is easy to verify that the first sum from (11) vanishes when P 2 = m 2 A) k 2 n = fi 2 . A hadron 
with the momentum k n must be off-mass-shell (the opposite case of an off-mass-shell recoil 
nucleus is irrelevant for the process under consideration). Due to the uncertainty principle, 
a hadron with the off-mass-shell momentum k n can be considered as a free particle during 
the time A ~ l/(fc n ,o — (^« + Z^ 2 ) 1 ^ 2 )- For fast particles, A is their propagation length. In our 
case A is similar for hadrons carrying k n and k. Thus, for virtual projectiles the contribution 
of the elastic channel, represented by the first sum from (11), is non- vanishing during the 
time of the order of A. If L is the duration of the scattering process, then the elastic cross 
section vanishes in the case (9), in agreement with the above euristic argumentation. In this 
case the only contribution to Im[f(s,t = 0)] is due to inelastic channel because the states 
\i > contain new particles, which can be on the mass shell. When (9) is fulfilled, we obtain 

Im[f A (s,t = 0)] = saf nel , a^O. (12) 

In electroproduction on nuclei, A ~ and La = where q is the virtual photon 
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momentum. In the shadowing region La > 100A and we have to use o A nel instead of a A ot to 
calculate F A {x). In most cases, the condition (9) is also valid for separate nucleons and the 
free nucleon structure function determines of nel rather than af ot for virtual hadrons. Note 
that the condition for the Fraunhofer diffraction never holds in electroproduction. 

In the independent particle approximation and for A ^> 1, the Glauber result for the 
nuclear inelastic cross section is given by 

'db(l-e-<« T ®). (13) 

In contrast to (5), the last equation is correct for any value of Im[f ]/ Re[f ]. In the 
exponent in (13) there is no 1/2-factor, which is present in (5) due to diffraction effects. If 
X/ta is finite, for the square well form of nuclear density the equation (13) can be rewritten 

as 

vLd = ZnpvtLi d W dzexp[-paf nel zexp(-z/X)}, (14) 
Jo Jo 

where p is the nuclear density. The exponent inside the exponent simulates the x-dependence 
of nuclear shadowing due to the finiteness of A/r^. To calculate nuclear structure functions 
F 2 ' VMD (x) at low x, we have used (3) with a to t equal to o A nel from (14). 

Following the data for high energy hadron-nucleon interactions, we assume that for virtual 
mesons o^ nel ~ 5/ Qa^ t as well as for free mesons. The final results are not sensitive to small 
variations of cr^ nel . We have used the standard set of parameters for low lying mesons. 
II(/i 2 ) for the meson continuum was determined from the data for the e + e~-annihilation 
following the prescription of Ref . [ D| . The cross section for the continuum was taken in 
a form <jf nel = 18mb GeV 2 / p 2 in order to reproduce the experimental nucleon structure 
function. For every x-bin, the average Q 2 is the same as in Ref.[[l[]. At x > 0.175, when the 
contribution of virtual hadrons becomes negligible, we have calculated F^' EMC (x) using the 
results of Ref. || with e=35 MeV, with the Fermi momentum p^=230 MeV/c and with the 



flux factors |19[ included. Our result for the structure function ratio can be represented as 

R = {F A {x) + AF 2 A ^(x))/AF 2 N (x), (15) 

where -F 2 A (x) = 0(0.175 — x)F A ' GVMD (x) + 9{x — 0.175)F^ 4,£;MC '(a;). The excess pion con- 
tribution AF A ' n (x) was calculated following the preceeding section with e= 35 MeV and 



with the pion three- momentum distribution from Ref. [12 . In some previous calculations 



on-mass-shell pions with E n = Jp%. + m\ || or static pions with E n = [[J have been used. 



Here we take the intermediate value E w = yp^ + ml/2. In this case R deviates by less than 
2% from the ratio calculated for the two extreme values of E n . For the number of excess 
pions per nucleon we obtained n n ~ 0.07. The elementary structure functions were taken 
from Ref.|2{J. 

The resulting R is shown by the solid line on Fig.l. The nucleon contribution, represented 
by the first term from (15), is shown by the dashed line. The earlier GVMD model result fl3[ 



for the nuclear shadowing effect is shown by the dotted line. In our method the nuclear 
shadowing is significantly larger than in previous calculations, because the shadowing is more 
effective in the inelastic channel. The agreement with the data is restored by the virtual 
pion contribution. In this case the antishadowing at x ~ 0.1 is also explained. Though 
we obtained a good description of the data, we understand that the GVMD model may be 
too simplified and can hardly provide a correct Q 2 -dependence of nuclear shadowing. But 
the result of this section is more general then the GVMD model: following our conjecture, 
the nuclear shadowing will be larger in other hadron dominance models based on Glauber 
theory. The agreement with data, obtained in some papers |T3|, |T4]| , can be recovered if the 
excess pion contribution will be included. 

4. The DY production by protons probes the antiquark content of a target and for the 
cross sections of this reaction it was usually accepted 

R DY = a A Y /Aa N DY » q A {x)/Aq N {x). (16) 

The measured i?Dy[0 is close to 1 (see Fig. 2). This would mean that there is no excess 
pions if (16) is correct. In Ref.|^] it is shown that there is no reason to neglect absorptive 
initial state (IS) interactions of fast partons in nuclei. The conclusion of that paper is that 
the color transparency forecasts restrict the strength of final state interactions but do not 
rule out IS interactions in hadropro duct ion. With the absorptive IS interactions included, 
for Rdy we obtain^] 

R DY = Eq A (x)/Aq N (x). (17) 



The factor F describes the projectile parton flux attenuation in nuclei, 

F = -± r [d 2 b(l-e-< b * T ®), (18) 

where a q bs is the cross section of absorptive interactions of fast initial partons with nucleons. 
We assume that a q abs is the same for quarks and antiquarks. In this case the parameter 
a q bs can be found from the DY productions by pions, which is not affected by the excess 
pion contribution. The analysis of Ref.[|J gives a q abs ~ 2mb and we used this value to 
calculate Rdy with (17) and (18). In our calculations, the nuclear antiquark density q A (x) 
includes the excess pion contribution by analogy to (1). This contribution was calculated 
as described in preceeding sections. We considered only u-quarks and antiquarks with the 



parametrizations from Ref.|20[. The binding correction for e=35 MeV was introduced as 
described above. The ratio Rdy calculated without the excess pion contribution (the lower 
dashed line on Fig. 2) is significantly below the data points. The full result for Rdy calculated 



with E-n = \jp^ + m ^/2 is shown by the solid line. The same for E v = is shown by the 
upper dashed line. As in the case of electroproduction, the pionic contribution improves the 
agreement with data. The bulk of the pion contribution is fixed by the sum rule (2) but the 
detailed form of resulting curves depends on the choice of the pion distribution. 

We conclude that both electro- and hadroproduction data, as well as the data for K + — A 



scattering can be explained with the excess pion contribution included. 
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Figure captions 

Fig.l. The ratio F^(x)/AF^(x) calculated without (dashed line) and with (solid line) 
excess pions. The dotted line is the earlier result of the GVMD model [0. The data are 
from Ref.0. 

Fig. 2. The ratio <x^ y / 'Aa^y for the Drell-Yan production by protons calculated without 
excess pions (lower dashed line) and with two different forms of excess pion distribution 
(solid and upper dashed lines). The earlier estimation of this ratio with excess pions included 
(dotted line) and the data points are from Ref. [|J. 
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